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We report recent results for lepton flavor violating (LFV) r decays obtained at Belle and 
BaBar with the world-largest data samples. 



1 Introduction 

An observation of the lepton flavor violation (LFV) is a clear signature of New Physics (NP) 
since the LFV in charged leptons has a negligibly small probability in the Standard Model (SM) 
even taking into account neutrino oscillations. For example, the branching fraction for the r 
decay into [i and 7, should be proportional to squared mass difference between neutrinos and 
the weak boson mass squared and is calculated to be less than 1 o -54--40 .[H Since the r lepton 
can have many possible LFV decay modes and is expected to be sensitive to NP because it has 
the heaviest mass among charged leptons, we search for the LFV from the r lepton as a NP 
phenomenon. 

The model including supersymmetry (SUSY), that is the most popular as the one beyond the 
SM, can naturally induce LFV at one loop. In many SUSY models, a r — > decay is expected 
to have the largest branching fraction among all the r LFV decays while a branching fraction 
for a r decay into \ir\ or can become the largest in some cases, such as the higgs-mediated 
case. By evaluating the branching fractions for various r LFV decays, we can distinguish what 
model is favored. Here, we especially focus on the three major modes, r — > ^7, ££'£" and £P° 
where £ = e or [i and P° = 7r°, r\ or rj ' . 

The world-highest r data samples has been collected in two B-factory experiments, by the 
Belle and BaBar collaborations, because the cross section for bb is very close to that for r + r~ 
on the T(45) mass. Therefore, the B-factory is a good place to search for the r LFV. Here, 
recent results with 0(1O 9 ) r samples obtained by Belle and BaBar are reported. 

2 Analysis Method 

In the r LFV analysis, in order to evaluate the number of signal events, two independent variables 
are defined, that are signal-reconstructed mass and energy in the center-of-mass (CM) frame 
from energies and momenta for the signal r daughters. In the r — > ^7 case, they are defined as 
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where En~ (P^) is a sum of the energies (a magnitude of a vector sum of the momenta) for 
\x and 7, the superscript CM indicates that the variable is defined in the CM frame and the 
E'beam means the initial beam energy in the CM frame. Principally, M^ 7 and AE should be m T 
(~ 1.78 GeV/c 2 ) and (GeV), respectively, for signal events while M^ 7 and AE will smoothly 
vary without any special structure in the background (BG) events. (BaBar often takes an 
energy-constrained mass (Mec) rather than the usually reconstructed mass.) The signal MC 
distribution of r — > fij on the M^ 7 - AE plane is shown in Fig. [TJ Due to the resolution, the 




Figure 1: (a) Distribution for the r — > signal events, (b) Comparison for the resolutions of the reconstructed 
t mass for r — s* /17 (blue), r — > fir/ (red) and r — > (black). 



signal events are distributed around M M7 ~ m T and AE ~ (GeV). Taking into account the 
resolution, we set the elliptic signal region. (See Fig. [TJ) Finally, we evaluate the number of 
signal events in the signal region. Generally, when some gammas are included, the resolution 
become worse. Therefore, r — > ££'£" has a good resolution while r — > £~f has a bad resolution. 
This means that good (bad) S/N is expected in the r — > ££'£" (r — > £7) analysis because the 
number of the BG events can be roughly estimated to be proportional to the size of the signal 
region. 

To avoid any bias for our analyses, we perform the blind analysis: Before fixing the selection 
criteria and the evaluation for the systematic uncertainties, we cover the data events in the 
signal region. 



3 r £-y 



In 2008, Belle has updated their result for r — > £7 with 4.9 x 10 8 r + r _ -pairs. ^1 The main 
BG events come from r —> £vv + extra 7 events and the radiative di-muon or bhabha events. 
The number of the signal events is evaluated by a generalized extended unbinned maximum- 
likelihood to the Mn^-AE distribution. As a result, the upper limit on the branching fraction 
for r — > /X7 (e-y) is set to be 4.5 x 10~ 8 (1.2 x 10~ 7 ) at 90% confidence level (CL). Recently, BaBar 
has published the new result with 4.8 x 10 8 r + r _ -pairs data sample on T(4S') as well as T(2S') 
and T(3»S) masses. New requirements to especially reject radiative di-muon or bhabha events 
are introduced. In the 2a signal region, The signal detection efficiency for r — > /X7 (e^f) is 6.1% 
(3.9%) while the expected number of the BG events is 6.1 ±0.5 (3.9±0.3) events. Consequently, 
2 (0) events are found in the r —> fij (eq) analysis and the 90% upper limit on the branching 
fraction is evaluated to be 4.4 (3.3) x 10~ 8 . Here, the 90% CL upper limits are calculated by a 
counting method. 
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Figure 2: Resulting distributions for t — !> ///y of Belle (a) and BaBar (b). The black dots and the shaded boxes 
show the data and signal MC distributions, respectively, and the (most inside) elliptic region is the signal region. 

4 t — y It' I" 

In the last summer, Belle and BaBar have shown the results for r — > ££'£" with their almost full 
data samples, of 7.2 x 10 s r + T~-pairs by Belle^and 4.3 x 10 8 ones by BaBar ^ . BaBar improves 
the efficiency for the lepton-ID and the better rejection for the background is obtained. As a 
result, no event in the signal region is found for all modes and the 90% CL upper limits on the 
branching fractions are set as shown in Table [TJ On the other hand, Belle takes very similar 
selection criteria to those for the previous analysis. Finally, for all modes, no event is observed 
in the signal region and the upper limits are evaluated around 2 or 3 times more sensitively than 
those obtained previously. 
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Table 1: Summary for the efficiency (Eff.), the expected number of the BG events, (iV^S) and the upper limit 
on the branching fraction (UL) for each mode. The upper and lower lines show the Belle and BaBar results, 

respectively. 
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Figure 3: Resulting distributions for r —¥ U'£" of Belle. The black dots and the shaded boxes show the data and 
signal MC distributions, respectively, and the elliptic region is the signal region. 



5 T — ^ £P° (P° = TT ,//,?/) 



In 2007, Belle and BaBar have published their results of the search for the r decay into a lepton 
and a neutral pseudoscalar (ir°, r/, rf). They have used around 400fb _1 and set the range of the 

|£?||M 

upper limits of (0.8 - 2.4) x 10~ 7 at 90% CL.^LII This summer, Belle reports the new result 
with about 900fb~ 1 . By studying the backgrounds in detail, they obtain about 1.5 times larger 
efficiency than previously in average and the expected number of the BG events in the signal 
region is achieved to be less than one for each mode. Details are summarized in Table. [2 Finally, 
one event is found in r — > ery(— > 77) while no event is observed in other modes. (See Fig. S|) 
The evaluated 90% upper limits on the branching fraction are (2.2 — 4.4) x 10 -8 . 



Mode (r ->■) 


Eff.(%) 


„oxp 
JV BG 


UL (xl0~ s ) 


Mode (r ->■) 


Eff.(%) 


„oxp 


UL (xl0~ s ) 


w(-> 77) 


8.2 


0.63 ±0.37 


3.6 


e^H 77) 


7.0 


0.66 ±0.38 


8.2 


/i?7(— > 7T7r7r u ) 


6.9 


0.23 ±0.23 


8.6 


e»7(— > 7r7T7r u ) 


6.3 


0.69 ±0.40 


8.1 


/i?7(comb.) 






2.3 


er^comb.) 






4.4 


[AT] (— > nitr)) 


8.1 


00 +U Ib 
u - uu -o.oo 


10.0 


er]'(— > irnri) 


7.3 


0.63 ±0.45 


9.4 


w'(-> 7P U ) 


6.2 


0.59 ±0.41 


6.6 


ev'(-> 7P°) 


7.5 


0.29 ±0.29 


6.8 


fir]' (comb.) 






3.8 


er;'(comb.) 






3.6 




4.2 


0.64 ±0.32 


2.7 




4.7 


0.89 ± 0.40 


2.2 



Table 2: Summary for the efficiency (Eff.), the expected number of the BG events, {N^q) and the upper limit on 
the branching fraction (UL) for each mode, where (comb.) means the combined result from subdecay modes. 
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Figure 4: Resulting distributions for r — > £P°. The black dots and the yellow boxes show the data and signal MC 
distributions, respectively, and the elliptic region is the signal region for each mode. This is preliminary. 



6 Summary 

We report recent results of the various search for LFV r decays obtained by Belle and BaBar. 
These results are about 100 times more sensitive than those by CLEO. This is achieved by not 
only much larger data samples but more effective BG rejection after detailed examination of the 
BG. At the present, we are updating results using full data samples accumulated by B-factories. 
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